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6 Abstract 

7 We determine the extinction-to-backscatter (S a ) ratios of dust using (1) airborne in-situ 

8 measurements of microphysical properties, (2) modeling studies, and (3) the Cloud- 

9 Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) observations 

10 recorded during the NASA African Monsoon Multidisciplinary Analyses (NAMMA) 

11 field experiment conducted from Sal, Cape Verde during Aug-Sept 2006. Using 

12 CALIPSO measurements of the attenuated backscatter of lofted Saharan dust layers, we 

13 apply the transmittance technique to estimate dust S a ratios at 532 nm and a 2-color 

14 method to determine the corresponding 1064 nm S a . This method yielded dust S a ratios of 

15 39.8 ± 1.4 sr and 51.8 ± 3.6 sr at 532 nm and 1064 nm, respectively. Secondly, S a at both 

16 wavelengths is independently calculated using size distributions measured aboard the 

17 NASA DC-8 and estimates of Saharan dust complex refractive indices applied in a T- 

18 Matrix scheme. We found S a ratios of 39.1 ± 3.5 sr and 50.0 ± 4 sr at 532 nm and 1064 

19 nm, respectively, using the T-Matrix calculations applied to measured size spectra. 

20 Finally, in situ measurements of the total scattering (550 nm) and absorption coefficients 

21 (532 nm) are used to generate an extinction profile that is used to constrain the CALIPSO 

22 532 nm extinction profile and thus generate a stratified 532 nm S a . This method yielded 

23 an S a ratio at 532 nm of 35.7 sr in the dust layer and 25 sr in the marine boundary layer 

24 consistent with a predominantly seasalt aerosol near the ocean surface. Combinatorial 

25 simulations using noisy size spectra and refractive indices were used to estimate the mean 

26 and uncertainty (one standard deviation) of these S a ratios. These simulations produced a 

27 mean (± uncertainty) of 39.4 (± 5.9) sr and 56.5 (± 16.5) sr at 532 nm and 1064 nm, 

28 respectively, corresponding to percent uncertainties of 15% and 29%. These results will 

29 provide a measurements-based estimate of the dust S a for use in backscatter lidar 

30 inversion algorithms such as CALIOP. 



31 


1. Introduction 


32 Lidar is a powerful tool for studying the vertical distribution of aerosols and clouds in the 

33 atmosphere. Of particular importance is the distribution and transport of Saharan dust 

34 systems. The deployment of CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder 

35 Satellite Observations), a joint NASA-CNES satellite mission, has enabled vertically 

36 resolved measurements of Sahara air layer(s) (SAL) which will provide significant 

37 insights into properties of Sahara dust aerosols. CALIPSO is designed to provide 

38 measurements to advance our understanding of the role of aerosols and clouds in the 

39 climate system [ Winker et al., 2009], The Cloud- Aerosol Lidar with Orthogonal 

40 Polarization [CALIOP, Winker et al., 2007] is the primary instrument on the CALIPSO 

41 satellite. CALIOP is designed to acquire vertical profiles of elastic backscatter at two 

42 wavelengths (1064 nm and 532 nm) from a near nadir- viewing geometry during both day 

43 and night phases of the orbit. In addition to the total backscatter at the two wavelengths, 

44 CALIOP also provides profiles of linear depolarization at 532 nm. Accurate aerosol and 

45 cloud heights and retrievals of extinction coefficient profiles are derived from the total 

46 backscatter measurements [Vaughan et al., 2009] The depolarization measurements 

47 enable the discrimination between ice clouds and water clouds [Hu et al., 2009] and the 

48 identification of non-spherical aerosol particles [Liu et al., 2009]. Additional information, 

49 such as estimates of particle size for the purpose of discriminating between clouds and 

50 aerosols, are obtained from the ratios of the signals obtained at the two wavelengths. On 

51 April 28, 2006, the CALIPSO satellite was launched into a low earth sun-synchronous 

52 orbit at a 705-km altitude, and an inclination of 98.2 degrees. A few months later, in 

53 August 2006, the NASA African Monsoon Multidisciplinary Analyses (NAMMA) 
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54 campaign commenced at the Cape Verde Islands, 350 miles off the coast of Senegal in 

55 West Africa. NAMMA was designed to study the evolution of precipitating convective 

56 systems, largely as this evolution pertained to the SAL and its role in the tropical 

57 cyclogenesis. Several aircraft flights were dedicated to nearly coincident measurements 

58 with NASA’s orbiting satellites (including Aqua, TRMM, and CloudSat/CALIPSO). For 

59 this study, we use data collected aboard NASA's DC-8 medium altitude research aircraft 

60 outfitted with, among other instruments, a full suite of sensors and probes designed to 

61 measure aerosol microphysical and optical properties. Relevant parameters include high 

62 spatial-resolution scattering and absorption coefficients at multiple wavelengths in the 

63 visible spectrum and dry particle size distributions over the 0.08 to 10 pm diameter range 

64 [ Chen et al. , 2010]. 

65 Depending on the mineralogical composition, the SAL can have a significant impact on 

66 both the radiation balance and cloud processes. Dust particles scatter in the shortwave 

67 regime (cooling the planet) and absorb both shortwave and longwave radiation (heating 

68 the planet). By some estimates the anthropogenic forcing due to dust is comparable to the 

69 forcing by all other anthropogenic aerosols combined [Sokol ik and Toon, 1996]. Saharan 

70 dust influences cyclone activity and convection in the region off the west coast of Africa 

71 and air quality as far west as the US east coast and Gulf of Mexico. There have been 

72 reports of causal links between cyclone activity and dust loading suggesting that perhaps 

73 the Sahara dust layer acts to inhibit cyclone development [Dunion and Velden, 2004] and 

74 more generally convection [Wong and Dessler, 2005]. Sahara dust is unique in its ability 

75 to maintain layer integrity as it is transported over long distances (-7500 km) to the 

76 Americas [Liu et al., 2008; Maring et al., 2003; Savoie and Prospero, 1976]. The 
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presence of Sahara dust layers have been found to perturb ice nuclei (IN) concentrations 

78 as far away as in Florida. During CRYSTAL-FACE (Cirrus Regional Study of Tropical 

79 Anvils and Cirrus Layers - Florida Area Cirrus Experiment), DeMott et al. [2003] found 

80 that IN concentrations were significantly enhanced in heterogeneous ice nucleation 

81 regimes warmer than -38°C, when Saharan dust layers are present. It is therefore 

82 important to study the distribution and optical properties of Sahara dust. 

83 In order to estimate the optical depth of the Sahara dust layers from elastic backscatter 

84 lidar measurements, the S a ratio must be known or prescribed. Given aerosol free regions 

85 above and below a lofted dust aerosol layer, S a can be calculated from the attenuated 

86 backscatter profile of a space-based lidar return [Young, 1995]. S a for dust aerosols is 

87 dependent on the mineral composition, size distribution, and shape parameters (e.g., 

88 aspect ratio, and complexity factor). All of these are highly variable and for the most part 

89 not well known. For these reasons, S a obtained from scattering models have larger 

90 uncertainties than models of the nearly spherical urban pollution or marine aerosols. 

91 There have been several studies and measurements of dust S a at 532 nm [. Ackermann , 

92 1998; Anderson et al., 2000; Berthier et al., 2006; Di Iorio et al., 2003; Di Iorio et al., 

93 2009; Muller et al., 2007; Muller et al., 2000; Tesche et al., 2009] and relatively few such 

94 measurements or studies of dust S a at 1064 nm [Ackermann, 1998; Liu et al., 2008; 

95 Tesche et al., 2009], Prior to NAMMA, there were a number of vertically-resolved 

96 measurements of Saharan dust microphysical and optical properties including AMMA 

97 and DODO campaigns NAMMA studies along with CALIPSO measurements provide a 

98 unique opportunity to compare extinction measurements derived from in situ profile 

99 measurements of total scattering and absorption aboard the NASA DC-8 and CALIPSO 
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100 extinction profiles estimated from two wavelength retrieval methods. These profiles by 

101 extension provide the constraints from which the lidar ratios can be determined as 

102 explained in the following sections. Section 2 discusses the CALIPSO lidar data and its 

103 analysis. The NAMMA data and analyses are discussed in Section 3, and coincident 

104 CALIPSO-NAMMA measurements are presented in Section 4. In Section 5, the size 

105 distributions measured during NAMMA aboard the DC-8 are implemented in a T-Matrix 

106 scheme to estimate profiles of S a ratios. Section 6 discusses the uncertainty in S a using a 

107 combinatorial method. 

108 

109 2. CALIPSO Lidar Data and Extinction-to-Backscatter Ratio Retrieval 

110 Methods 

111 The CALIPSO lidar data used for these studies are the version 2.01 lidar level 1 

112 attenuated backscatter returns at the 532 nm perpendicular and parallel channels, and 

113 1064 nm total attenuated backscatter. The volume depolarization ratio is determined from 

114 the perpendicular and parallel channels and used to identify dust aerosols [Liu et al., 

115 2009; Omar et al., 2009]. For NAMMA underflights of CALIPSO and near spatial 

116 coincidences where both missions observed dust layers of optical depths greater than 

117 about 0.3, we compare the extinction profiles from in situ measurements to CALIPSO 

118 profiles. In such cases, we calculate the extinction using S a that was determined using the 

119 transmittance method or an S a ratio constrained by the in-situ extinction profiles. In both 

120 cases, we use the 2-color methods to retrieve the 1064-nm S a , after determining the 532- 

121 nm S a . These two methods, transmittance and 2-color, are discussed below. 
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2.1. Transmittance Methods 


123 The transmittance method uses the following equation describing the relationship 

124 between optical depth and integrated attenuated backscatter, as in Platt [1973]: 

125 

126 y' = — (l-exp(-2riT)) (1) 

127 

128 Here y' is the integrated (from layer base to top) attenuated backscatter, 

129 

(•top 

130 Y '= [ Pa( r )T a 2 (r)dr (2) 

J base 

131 

132 t is optical depth, q is a multiple scattering parameter, T 2 = exp(-2qx) is the layer- 

133 effective two-way transmittance, and S a = c a /p a where p a is the aerosol backscatter 

134 coefficient and a a is the aerosol extinction coefficient. This ratio is assumed constant 

135 throughout a feature. Note that the quantities S a , y’, and x describe characteristics of an 

136 aerosol layer, i.e., they are associated with the backscatter and extinction of aerosol 

137 particles only. If we define an effective S a ratio, S* = qS a , we can rewrite Eq. (1) as 

138 follows: 

1-T 2 

139 S* = (3) 

2y' 

140 

141 The effective two-way transmittance is typically obtained by fitting the returns both 

142 above and below a feature to a reference clear air scattering profile obtained from local 
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143 rawinsonde measurements or meteorological model data [Young, 1995], In this study, the 

144 transmittance method is used to determine S* from the 532-nm CALIPSO measurements 

145 whenever clear air scattering signals are available both above and below an aerosol layer. 

146 However, the same method is not applicable to the 1064-nm CALIPSO measurements, 

147 because a reliable measurement of the clear air scattering at 1064 nm, which is about 16 

148 times smaller than that at 532 nm, is not available. To determine S* at 1064 nm, the 2- 

149 color method described in the next subsection is used. 

150 

1 5 1 2.2. The 2-Color Method 

152 The 2-color or two- wavelength method was first proposed by Sasano and Browell 

153 [ Sasano and Browell, 1989] and adapted to space-borne lidar measurements using an 

154 optimization technique by Vaughan et al. [2004], The method requires apriori knowledge 

155 of S a at 532 nm and a suitable profile of 532-nm attenuated backscatter amenable to the 

156 calculation of 532-nm aerosol backscatter coefficient profiles. For the NAMMA cases 

157 described below these preconditions were satisfied. Whenever a suitable region of clear 

158 air was identified both above and below an aerosol layer, the S a at 532 nm was 

159 determined using the transmittance method described above. Here clear air layer is 

160 defined as a region of low attenuated scattering ratios with a mean value equal to or less 

161 than 1 and a slope with respect to altitude of approximately zero. This is further 

162 confirmed by low volume depolarization ratios in a small region (~ 1/2 km) below the 

163 aerosol layer. In cases where coincident NAMMA measurements are available, S a ratio at 

164 532 nm is the value that provides the best fit between the retrieved CALIPSO extinction 

165 profiles and the NAMMA in situ extinction profiles obtained by summing the total 
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scattering and absorption measured by a nephelometer and a Particle Soot/ Absorption 
Photometer (PSAP), respectively, aboard the DC-8. 


167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 


Once S a ratio is determined at 532 nm, the value at 1064 nm can be calculated using the 
2-color method. We note that this technique can be used to derive S a at 532 nm if the 
value at 1064 nm is known. Given a solution of the particulate backscatter at 532 nm, 
P 532 , p , the 2-color method uses a least squares method to minimize the difference between 
the measured attenuated total backscatter at 1064 nm, Bio 64 , and the attenuated 
backscatter at 1064 nm (right hand side of eq. (4)) reconstructed from the extinction and 
backscatter coefficients at 532 nm. 


®10<S4 ( r ) ( Pm, 1064 ( r ) + Pp,1064 ) ' ^p,1064 ( f ) 

= (P m,1064 (0 + X ‘ Pp432 (*)) ' (~2 ' • X ' Y 532 (f)) 


( 4 ) 


The only unknowns (underlined) in eq.(4), are the S a ratio at 1064 nm, Sio 64 , and the 
backscatter color ratio, % (defined as Pio 64 , p /P 532 , p )- These are both intensive aerosol 
properties defined by the layer composition, size distribution, and shape of its constituent 
particles. Since these characteristics do not vary substantially in a given aerosol layer, we 
make the assumption that Sio 64 and % are constant within the layer. The algorithm details 
and optimization techniques are discussed at length in Vaughan (2004) and Vaughan et 
al. [2004] 
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186 3. Numerical Calculation Based on NAMMA In-situ Measurements 

187 3.1. Aerosol Microphysical Properties 

188 We use measurements of the aerosol size distributions based on number from the 

189 Aerodynamic Particle Sizer (APS, TSI Incorporated, Shoreview, MN) and the Ultra-High 

190 Sensitivity Aerosol Spectrometer (UHSAS, Droplet Measuring Systems, Boulder, CO). 

191 The UHSAS measures the fine mode aerosol size distributions from 0.06 to 0.98 pm, and 

192 the APS measures the coarse mode size distributions from 0.6 to 5.5 pm. We use the size 

193 distributions to identify the presence of aerosol dust layers. In many cases, these intense 

194 dust layers were visually identified by the instrument operators and in some cases 

195 specifically targeted by the DC-8 operators for sampling. Additional information about 

196 the composition of these layers is available from the NAMMA data archives 

197 (http://namma.msfc.nasa.gov/). For each size distribution sampled during a 5-second 

198 interval we fit the discrete measurements to the best continuous bimodal lognormal size 

199 distribution, as shown in Figure 1. The geometric mean radius and standard deviation of a 

200 fine and coarse mode derived from the in situ measurements are used in the numerical 

201 calculations. 

202 

203 On August 25, 2008, the DC-8 flew through a dense elevated dust layer measuring nearly 

204 1 km in thickness at a mean altitude of about 2.3 km in an hour-long mostly straight and 

205 level flight. The APS and UHSAS measured coarse and fine size distributions, 

206 respectively, through the dust layer at intervals of 5 seconds. Figure 2 shows the 

207 probability distributions of the geometric mean fine and coarse radii of the dust layer. 

208 This figure is generated by taking the discrete 5-second size distribution measurements 
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209 and fitting these to a bimodal lognormal distribution as described above. For this dust 

210 layer, the microphysical properties of mean, median, and standard deviations of the fine 

211 ( coarse ) radius distributions, as shown in Figures 2 (a) and (b), are 0.059, 0.061, and 

212 0.0064 pm {0.54, 0.57, and 0.083 jum), respectively. The mean, median, and standard 

213 deviations shown in Figure 2 (c) and (d) for the fine {coarse) geometric standard 

214 deviation distributions are 1.613, 1.630, and 0.101 {1.495, 1.545, and 0.151), 

215 respectively. 

216 

217 The distributions in Figure 2 show that dust properties after lofting of these layers remain 

218 relatively unchanged. Other studies [Liu et al., 2008; Maring et al., 2003; Prospero and 

219 Carlson, 1971; 1972], have shown the same consistency in properties after long range 

220 transport of dust. In each case, the means and medians of the size distribution descriptors 

221 are close, i.e., the size descriptors are nearly normally distributed. The standard deviation 

222 is a small fraction (<0.16) of the means, i.e., the variance of the data is small and thus the 

223 layer is quite homogenous with respect to size across the 1 km vertical extent of the dust 

224 plume. 

225 

226 3.2. Scattering Models 

227 Mie scattering calculations [ Mie , 1908], when applied to dust, are adequate for total 

228 scattering, albedo and other flux related quantities, but result in large errors when used to 

229 retrieve optical depth from satellite reflectance measurements. In particular and central to 

230 the theme of this paper, are S a ratios calculated from measured size distributions. Mie 

231 calculations underestimate S a by up to a factor of 2.0 leading to substantial errors in the 
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lidar derived aerosol optical depths [Kalashnikova and Sokolik, 2002], This has been 

233 known experimentally for quite some time: laboratory measurements by Perry et al 

234 [1978] showed non-spherical particles, when compared to spherical particles having the 

235 same equivalent volume, enhance side scattering and suppress backscattering. To account 

236 for non-sphericity of dust particles, we use T-matrix calculations with the assumption that 

237 the dust shapes can be modeled by randomly oriented prolate spheroids. T-Matrix is a 

238 matrix formulation of electromagnetic scattering first proposed by Waterman [1971], and 

239 subsequently improved and extended to much larger sizes and aspect ratios by Mischenko 

240 et al. in a series of papers [Mishchenko, 1991; 1993; Mishchenko and Travis, 1994; 

241 Mishchenko et al., 1996a; Mishchenko et al., 1996b; Mishchenko and Travis, 1998]. The 

242 T-Matrix code used in these calculations is described in detail in Mishchenko and Travis 

243 [1998], This method is particularly suitable for light scattering calculations of non- 

244 spherical, polydisperse, randomly oriented particles of identical axially symmetric shape 

245 with size parameter, x (x= 7 id p /X; d p is particle diameter and X is the wavelength), smaller 

246 than 30. 

247 

248 It is challenging to determine representative statistics of the mean shape for dust particles 

249 because of their complexity and variety in shape. These particles are not only confined to 

250 desert regions but are ubiquitous in continental areas where they contribute quite 

251 significantly to the extinction budget [Omar et al., 1999], Fortunately, however, 

252 randomly oriented prolate and oblate spheroids can adequately represent the scattering 

253 properties of non-spherical particles of the same aspect ratio [cf. Mischenko et al., 1995; 

254 Mishchenko and Travis, 1994], The aspect ratio is the ratio of the largest to the smallest 
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255 particle dimension. A prolate (oblate) spheroid is a rotationally symmetric ellipsoid with 

256 a polar diameter greater (smaller) than the equatorial diameter. 

257 

258 There have been several studies of aspect ratio distributions representative of dust 

259 aerosols. From an analysis of scanning electron microscope images of yellow sand 

260 particles, Nakajima et al. [1989] found that the distribution of the minor to major particle 

261 radius ratio peaked around 0.6, equivalent to an aspect ratio of 1.67. An investigation of 

262 mineral dust particle shapes using electron microscopy by Okada et al. [1987] found a 

263 mean aspect ratio of 1.4 ranging from 1.0 to 2.3. Hill et al [1984] compared the measured 

264 scattering properties of 312 samples of soil dust with the simulated scattering properties 

265 of randomly oriented prolate spheroids using T-matrix. They found the distributions of 

266 the aspect ratio of prolate spheroids centered near 2.3 most closely reproduced the 

267 measured scattering properties. For this study, we use a mean aspect ratio of 2.0 based on 

268 the above studies and investigate the sensitivity of the S a ratio to aspect ratios ranging 

269 from 1.7 to 2.3 partly to account for the dependence of the aspect ratio on size as 

270 discussed by Kalashnikova and Sokolik [2004], 

271 

272 There is quite a wide range of estimated and measured mineral dust refractive indices (m- 

273 ik). For wavelengths of 550 and 1000 nm, d’Almeida et al. [1991] estimate the real part 

274 (m) of 1.53 and 1.52, respectively, and a spectrally invariant imaginary part (k) of 0.008 

275 for dust-like aerosols, and 1.53-i0.0055 and 1.53-i0.001, respectively, for mineral dust. 

276 Ackerman (1998) used values of 1.53-i0.0043 and 1.53-i0.0063 to calculate dust S a ratios 

277 of 19-23 sr and 17-18 sr at 532 and 1064 nm, respectively. These values are much lower 


12 



278 than more recent 532-nm dust S a ratios of 40-60 sr [Cattrall et al., 2005; Di lorio et al., 

279 2003; Muller et al., 2007; Murayama et al., 2003; Voss et al., 2001] because the 

280 Ackerman study assumed spherical particles. 

281 Retrievals from radiances measured by ground-based Sun-sky scanning radiometers of 

282 the Aerosol Robotic Network (AERONET) over a 2-year period yielded dust complex 

283 refractive index values of 1.55 ± 0.03 - i0.0014 ± 0.001 at 670 nm and 1.55 ± 0.03-i0.001 

284 ± 0.001 at 1020 nm at Bahrain (Persian Gulf), and 1.56 ± 0.03 - i0.0013 ± 0.001 at 670 

285 nm and 1.56 ± 0.03-i0.001 ± 0.001 at 1020 nm at Solar Village in Saudi Arabia [Dubovik 

286 et al., 2002], Using vertically resolved aerosol size distributions in a scattering model 

287 constrained by lidar measurements of aerosol backscattering coefficient at 532 nm, Di 

288 lorio et al. (2003) estimated dust refractive indices of 1.52 to 1.58 (real part), and 0.005 

289 to 0.007 (imaginary part). Kalashnikova and Sokolik [2004] calculated effective 

290 refractive indices from component mixtures and found values of 1. 6140.0213 and 1.59- 

291 i0.0032 for Saharan dust at wavelengths of 550 nm and 860 nm, respectively. The values 

292 for Asian dust, from the same monograph, are 1.5140.0021 and 1.5 140. 0007 at 550 nm 

293 and 860 nm, respectively. Using a twin angle optical counter , Eidehammer et al.(2008) 

294 estimated the indices of refraction to be in the range 1 .60-1 .67 for the real part and 0.009- 

295 0.0 1 04 for the imaginary part. 

296 

297 The Saharan Mineral Dust Experiment [SAMUM, Heintzenberg, 2009; Rodhe, 2009] 

298 based in Morocco in 2006, produced several independent estimates of the complex 

299 refractive indices of Saharan dust. Kandler et al. [2009] determined Saharan dust aerosol 

300 complex refractive index from chemical/mineralogical composition of 1.5540.0028 and 
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301 1.57-i0.0037 at 530 nm for small (diameter < 500 nm) and large particles (diameter > 500 

302 nm), respectively. Schladitz et al. [2009] derived mean refractive indices of 1.53 - 

303 i0.0041 at 537 nm and 1.53 - i0.0031at 637 nm from measurements of scattering and 

304 absorption coefficients, and particle size distributions. Using similar methods during 

305 SAMUM, Petzold et al. [2009] found real parts of the refractive indices of Saharan dust 

306 ranging from 1.55 to 1.56 and imaginary parts ranging from 0.0003 to 0.0052. Some of 

307 the estimates of refractive indices reported in the literature are summarized in Table 1 . 

308 
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309 Table 1. Summary of complex dust refractive indices from previous studies 


Wavelength 

(nm) 

Real Part 

Imaginary Part 

Source 

500 

1.50 

0.0045 

Volz [1973] 

1000 

1.53 

0.008 

d 'Almeida et al. [1991], dust- like 

550 

1.53 

0.0055 

d’ Almeida et al. [1991], mineral dust 

1000 

1.53 

0.001 

d’ Almeida et al. [1991], mineral dust 

532 

1.53 

0.0043 

Ackerman [1998] dust 

1064 

1.53 

0.0063 

Ackerman [1998] dust 

670 

1.55 ±0.03 

0.0014 ±0.001 

Dubovik et al. [2002], Bahrain dust 

1020 

1.55 ±0.03 

0.003 ±0.001 

Dubovik et al. [2002], Bahrain dust 

670 

1.56 ±0.03 

0.0013 ±0.001 

Dubovik et al. [2002], Solar Village dust 

1020 

1.56 ±0.03 

0.001 ±0.001 

Dubovik et al. [2002], Solar Village dust 

532 

1.52-1.58 

0.005-0.007 

Di Iorio et al. [2003], Saharan dust 

500 

1.42 

0.003 

Israelevich [2003], Sede Boker dust 

860 

1.51 

0.0032 

Kalashnikova and Sokolik [2004] Saharan dust 

550 

1.51 

0.002 

Kalashnikova and Sokolik [2004] Asian dust 

860 

1.51 

0.0007 

Kalashnikova and Sokolik [2004] Asian dust 

670 

1.45 

0.0036 

Omar et al. [2005] dust cluster 

550 

1.53 

0.0015 

Cattrall et al. [2005] mineral dust 

1020 

1.53 

0.0005 

Cattrall et al. [2005] mineral dust 

300-700 

1.60-1.67 

0.009-0.0104 

Eidehammer et al. [2008] Wyoming dust 

537 

1.53 

0.0041 

Schladitz et al. [2009] Saharan dust 

637 

1.53 

0.0031 

Schladitz et al. [2009] Saharan dust 

450 

1.55-1.56 

0.0003-0.0052 

Petzold et al. [2009] Saharan dust 

700 

1.55-1.56 

0.0003-0.0025 

Petzold et al. [2009] Saharan dust 

530 

1.55 

0.0028-0.0037 

Kandler et al. [2009] Saharan dust 


310 


311 To perform the sensitivity study described in section 6, we use values of the real part of 

312 the refractive index ranging from 1.45 - 1.55 (normally distributed) and the imaginary 
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313 part ranging from 0.00067 - 0.006 (log normally distributed) with a central value of 1 .50- 

314 i0.002. For the scattering calculations using NAMMA size distributions, we use the 

315 central values of the refractive indices along with the nearly instantaneous (5 second 

316 interval) size distribution measurements to generate profiles of the aerosol properties. 

317 Figure 3 is a plot of the fine mode, coarse mode, and total phase function of the dust 

318 plume encountered on August 19, 2006. The fine mode and coarse mode phase functions 

319 are computed from the mean of the instantaneous size distributions, and the total phase 

320 function is the area-weighted composite of the fine and coarse mode phase functions. The 

321 phase functions are driven largely by the coarse mode, especially at 532 nm, and exhibit a 

322 more pronounced peak in the backscattering direction at 532 nm than 1064 nm 

323 

324 4. Data Analyses of CALIPSO NAMMA Coincident Measurements 

325 For this study we analyzed coincident measurements of the CALIOP 532-nm extinction 

326 profiles and in situ extinction profiles measured at wavelengths near the CALIOP green 

327 channel. The in situ extinction coefficient is obtained by summing the scattering (550 

328 nm) and absorption (532 nm) coefficients measured by a nephelometer and a PSAP, 

329 respectively. Data from the nephelometer and PSAP have been corrected for errors 

330 associated with the limited detector viewing angle [. Anderson and Ogren, 1998] and 

331 scattering from the filter media [ Virkkula et al., 2005], respectively. We make the 

332 assumption that the scattering properties are invariant over the 532 - 550 nm range for 

333 these large dust particle sizes. We chose three days on which there were near collocated 

334 CALIPSO and NAMMA measurements of nearly the same airmass. 
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335 4.1. August 19, 2006 NAMMA Flight 4 

336 August 19 was one of the days the DC-8 performed an under flight of CALIPSO. Figure 

337 4 shows the time-altitude flight track of the DC-8. A nearly coincident in-situ profile was 

338 obtained during the second ascent leg shown in the figure. Atmospheric context for this 

339 flight is given by Figure 5(a) and (b), in which both the DC-8 and CALIPSO flight tracks 

340 are superimposed on images of measurements made by the Moderate Resolution Imaging 

341 Spectroradiometer (MODIS) and Measurements of Pollution in the Troposphere 

342 (MOPITT), respectively. The DC-8 flight tracks are the black irregular octagons and the 

343 CALIPSO orbit tracks are the straight lines in both images. The underflights in the 

344 figures are segments where the DC-8 flight tracks are nearly exactly collocated and 

345 parallel to the CALIPSO flight tracks. The MODIS optical depth near the coincident 

346 flight segment is about 0.5. The MOPPITT image shows moderate CO concentrations in 

347 the vicinity of the coincident flight track and therefore indicates that most of the aerosol 

348 is dust and not continental pollution or biomass burning. This can also be confirmed by 

349 the CALIPSO depolarization measurements. The image also shows high CO 

350 concentrations (>2.5x10 molecules/cm ) to the south of the DC-8 flight tracks most 

351 likely due to biomass burning, and identified as such by the CALIPSO aerosol subtyping 

352 scheme illustrated later in this section. 

353 

354 Unfortunately, the direct underflight of CALIPSO by the DC-8 was a nearly level flight 

355 with no in situ profile information. In fact the DC-8 was at a high altitude near Flight 

356 Level 330 (~ 10 km) throughout the underflight and therefore did not encounter any 

357 significant aerosol layer at this altitude to sample. We used the ascending leg of the DC-8 
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358 flight which corresponds to the flight segment during ascent to the CALIPSO underflight 

359 portion denoted by blue dots in Fig 5a. The in situ profiles segment is shown in Fig. 4. 

360 The CALIPSO browse images (e.g., Fig. 6) are plots of the attenuated backscatter color 

361 coded by intensity varying from blue (weak) to white (very strong). A horizontal line 

362 near the 0 km mark denotes the surface. The CALIPSO data used for comparison with the 

363 in-situ data come from the 80 profiles (~27 km horizontal average) between the white 

364 lines in the browse image (Fig. 6). To determine the optimal 532-nm S a ratio for these 

365 data, we iteratively adjust S a until the difference, in a least squares sense, between the 

366 retrieved CALIPSO profile, a a , 532 -nm, and the measured NAMMA profile, c a , 55 o-nm, is 

367 minimized. 

368 

369 Figure 7 is a plot of the extinction profiles retrieved from CALIPSO’s 532-nm 

370 backscatter profiles and the profiles of the sum of scattering (550 nm) and absorption 

371 (532 nm) measured by the nephelometer and the PSAP aboard the NASA DC-8, 

372 respectively. These profiles were taken during the ascent leg shown in Figure 4 

373 corresponding to the flight tracks shown in Figure 5(a). The boundary between the two 

374 layers is determined by the increase in the extinction coefficient near 2 km. The two S a 

375 ratios are values that provide the best fit in a least squares sense of the CALIPSO data to 

376 the NAMMA in-situ measurements. The 532-nm S a ratios are consistent with a dust 

377 plume (S a = 35.7 sr) above seasalt in the marine boundary layer (S a = 25 sr). The root 

378 mean square (rms) of the differences between the CALIPSO and the NAMMA in-situ 

379 extinction coefficient profiles are 20.2 Mm' 1 , and 47.4 Mm' 1 , for the marine and dust 

380 layers, respectively. The rms of the differences in the extinction coefficients in the clear 
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381 region (4 to 8.5 km) above the dust layer in Fig. 6 is 12.3 Mm' 1 . We see some differences 

382 in the two extinction profdes in Figure 7 likely due to the temporal and spatial differences 

383 of 30 minutes and 160 km, respectively, between DC-8 and CALIPSO. 

384 

385 Figure 8 is an image of the aerosol scattering ratios measured by the Lidar Atmospheric 

386 Sensing Experiment (LASE) on board the DC-8 during the ascent leg of Flight 4 on 

387 August 19, 2006 as indicated in Figure 4. The measurements were made at 815 nm and 

388 the extinction calculation was performed using an S a ratio of 36 sr [Ismail et al., 2010]. 

389 The DC-8 flight track is shown by the solid line in Figure 8 during which the in-situ 

390 measurements shown in Figure 7 were made. The figure also shows the region nearest to 

391 the CALIPSO overpass where measurements of the CALIPSO profdes shown in Figure 7 

392 were made. LASE also observed a dust layer extending to an altitude of 6 km at the 

393 coincident point. This dust layer is optically and geometrically thick, and is lofted over a 

394 layer of lower optical depth (aerosol scattering ratio ~ 1) in the marine boundary layer. 

395 Some differences in altitude and extent of the layer between the CALIPSO and LASE 

396 measurements can be attributed to temporal and spatial mismatch. 

397 

398 Figure 9 is a plot of the results of (a) the cloud-aerosol discrimination and (b) aerosol 

399 classification algorithms applied to the data shown in the browse image including the 

400 NAMMA underflight (Flight 4 of August 19, 2006). The CALIPSO level 2 algorithms 

401 first discriminate between aerosol and clouds [Liu et al., 2009] and then classify the 

402 aerosol layers into aerosol subtypes [ Omar et al., 2009], Figure 9(a) shows that some of 

403 the optically thick aerosol near 15° N was misclassified as clouds, and thus was not 
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404 examined by the aerosol subtyping algorithm. The presence of biomass burning smoke 

405 and polluted dust (mixture of dust and smoke) during the first part of the flight depicted 

406 in Figure 9(b) is borne out by the high CO concentrations in the MOPPITT data to the 

407 south of the DC-8 flight tracks shown in Figure 5(b). Though the small lump of aerosol at 

408 the surface near 1 5° N is classified as pure dust, it is more likely a mixture of dust and 

409 marine aerosol. The white line in the Figure is the midpoint of the 80 profiles averaged 

410 for the retrieval of the extinction profile discussed above. 

411 

412 CALIPSO also measured a dense Saharan dust layer to the southwest of the coincident 

413 measurements during a nighttime orbit on the same day (August 19, 2006). The browse 

414 images of attenuated backscatter at 532 nm for this measurement are shown in Figure 10. 

415 The inset map in Figure 10 shows the CALIPSO ground track in blue. This dust layer 

416 appears to be a more robust part of the same dust plume observed during the coincident 

417 measurement. 

418 

419 As shown (by the red dotted oval) in Figure 10(a), the layer exceeds 1000 km in 

420 horizontal extent (from 18.3N to 10. 3N). The 532-nm aerosol optical depth (AOD) is 

421 greater than about 0.3 across most of the 1000-km orbital segment shown in Figure 10(a). 

422 Figure 10(b) is a magnified illustration of the region in 9(a) subtended by the yellow 

423 dotted line. We divided the layer into five segments and applied the transmittance method 

424 of section 2.1 to calculate a 532-nm S a and the 2-color method of section 2.2 to calculate 

425 a 1064-nm S a ratio. These values are shown in green (532 nm) and red (1064 nm) in 

426 Figure 10(b). Clear air regions above and below each dust layer were identified 
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427 


manually, by inspection of the profiles. Note that for this mesoscale layer the 532-nm S a 

428 ratios range from 38 to 41 sr with an average of 40.1 sr and the 1064-nm S a ratios range 

429 from 45.8 to 54.2 with an average of 50.9 sr. 

430 

431 Figure 10 shows that the Sahara dust layers once elevated are consistent both 

432 geometrically (the layer is confined to 3-5 km altitude band) and optically (S a variation at 

433 both wavelengths is small and the layer optical depth is near 0.3). The S a ratio is an 

434 intensive aerosol property that depends on the composition, size distributions, and 

435 particle shape of the aerosol and its consistency is an indication that these layers stay 

436 intact over very long distances. Other studies have shown the transport of relatively 

437 unmixed Saharan mineral dust to the south American rainforest [Ansmann et al., 2009; 

438 Graham et al., 2003] and western Atlantic Ocean [Formenti et al., 2003; Kaufman et al., 

439 2005], including the US eastern seaboard (cf. Liu et al. 2008). 

440 4.2. August 26, 2006 

441 The August 26 DC-8 flight included an underflight of CALIPSO. As is the case with the 

442 coincident CALIPSO-NAMMA measurements on August 19, the collocated 

443 measurements are at one level and lack in-situ profile measurements. In situ profiles of 

444 the size distributions were estimated from the DC-8 data obtained during the descent 

445 flight segment shown by a red dashed tilted oval in Figure 11 (a) and flown about two 

446 hours earlier than the CALIPSO underflight. The MOPITT CO levels (Figure 1 lb) for 

447 this period are slightly elevated. The extinction comparison shows that the layer observed 

448 by CALIPSO is more elevated than the one encountered by the DC-8. Figure 12 is a 

449 browse image of the CALIPSO 532-nm attenuated backscatter coefficients measured 
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450 during the orbital segment corresponding to Flight 8 of the DC-8. The CALIPSO data 

45 1 used for comparison with the in-situ extinction profiles were extracted from the region 

452 between the two white lines shown in the figure. 

453 

454 Figure 13 compares the CALIPSO extinction profile with a profile of the extinction 

455 derived from in situ measurements (i.e., the sum of the absorption and scattering 

456 coefficients). The maximum extinction coefficients (160 and 140 Mm' 1 by CALIPSO and 

457 the DC-8, respectively) are comparable, showing that the layer is intact after two hours. 

458 The altitudes of the maximum layer extinctions for the CALIPSO and NAMMA 

459 measurements are offset by at least one kilometer. The in situ measurements and 

460 CALIPSO observations are far removed from each other in this case (two hours and 1250 

461 km). The data shown in Figure 12 is north of the NAMMA flight segment. These 

462 differences are shown by the mismatch in layer heights shown in Figure 13. Because of 

463 the relatively large mismatch, the constraint method, which requires good coincidence, is 

464 not applicable to this case to derive S a . However, the transmittance technique can be 

465 applied using the CALIPSO measurement averaged over the region bounded by two 

466 white lines in Figure 12. Note that the dust layer overlies a streak of marine stratus just 

467 above the marine boundary layer. At 532 nm, the dust layer has an optical thickness of 

468 0.41 and an S a ratio of 38.2 sr. Though the Saharan layers can have horizontal extents of 

469 thousands of km, it is possible that the layer observed by CALIPSO is not the same as the 

470 one measured by the in situ instruments aboard the DC-8. Notwithstanding this 

471 possibility, the properties measured by the two methods provide independent 

472 characterization of the Saharan dust layer(s). 
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474 4.3. September 1, 2006 

475 Flight 10 of the DC-8 on September 1, 2006 neither directly underflew nor intercepted 

476 the CALIPSO orbit tracks. However, this flight sampled some of the same aerosol layers 

477 measured during a CALIPSO orbit track, as shown in the MODIS aerosol optical depth 

478 image in Figure 14(a). The region of interest is denoted by the red circle. The DC-8 flight 

479 segment most relevant and of closest proximity is the descending leg located in this 

480 region. The MOPPITT CO concentrations are low (1.5 - 2 x 10 18 molecules/cm 2 ) in the 

48 1 comparison region. It is therefore likely that most of the aerosol mass is Saharan dust. 

482 CALIPSO preceded the DC-8 by 14 hours. As in the previous flight of August 26, the 

483 observed layers by CALIPSO and the in situ measurements may be different. 

484 

485 The CALIPSO browse image (Figure 15) shows the layers of interest marked by two 

486 white lines and comprising of 80 profiles. The MOPITT image of this area shows no 

487 enhancement of CO in the sampling region subtended by a red circle in Figure 14(a). The 

488 transmittance and 2-color methods yielded S a ratios of 39.8 sr and 56 sr at 532 nm and 

489 1064 nm, respectively. The optical depth determined from the inversion of the CALIOP 

490 attenuated backscatter using this lidar ratio (39.8 sr at 532 nm) is 0.55. In Figure 14(a), 

491 the MODIS AOD is about 0.5 in this region. 

492 

493 It is a general transport pattern that Saharan dust aerosol goes through a phase of rising 

494 motion near the source then relative horizontal transport culminating in descent and 

495 deposition near the Americas [ Ansmann et ah, 2009; Formenti et al., 2003; Graham et 
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496 al., 2003; Kaufman et al., 2005; Okin et al., 2004], If the layer observed on September 1, 

497 2006 by CALIPSO is the same as the one observed by the in situ measurements, then it 

498 was in the rising phase at a rate of approximately 0.1 kilometers per hour. The rate of 

499 ascent is based on the time difference between the CALIPSO overpass (3:30 am Local 

500 Time) and DC-8 flight track (6:30 pm Local Time) of nearest approach. The extinction 

501 profiles’ comparison (Figure 16) between the DC-8 in-situ measurements and the 

502 CALIPSO measurements shows an offset of 1.5 km in the altitude of maximum 

503 extinction coefficient. The layer shown at the same location has been elevated by 1.5 km 

504 since CALIPSO sampled it. During this rising phase, there is no evidence of significant 

505 deposition, since the maximum extinction coefficient, a property of the aerosol loading, 

506 does not decay appreciably. 

507 

508 5. Extinction-to-backscatter (S a )ratios calculation based on NAMMA in situ 

509 measurements 

510 To determine profiles of S a ratios and validate the retrieved values, we perform numerical 

511 calculations for Saharan dust based on NAMMA in-situ size distribution measurements. 

512 We use the DC-8 APS and UHSAS measurements (Chen et al., 2010) to determine coarse 

513 and fine size distributions, and then calculate coarse and fine mode phase functions, as in 

514 Figure 3, using a T-Matrix scheme. The S a ratio of the aerosol is derived from an area- 

515 weighted integral of the two modes. Figure 17 is a plot of the profile of S a ratios of the 2- 

516 km dust layer encountered by NAMMA Flight 4 on August 19, 2006. The figure shows a 

517 profile of the coarse number concentration which marks the bottom and top of the layer at 

518 2.5 and 4.6 km, respectively. This is very similar to the coincident CALIPSO extinction 
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519 profile shown in Figure 7. The 532-nm and 1064-nm S a ratios calculated by this method 

520 are 34.3 ± 2.0 sr and 50.2 ± 5.7 sr, respectively. These are in good agreement with S a 

521 ratios of 38 to 41 sr at 532 nm and 45.8 to 54.2 sr at 1064 nm independently determined 

522 from CALIPSO measurements using the transmittance technique (Figure 10) on the same 

523 day, albeit further downfield. 

524 Part of the flight on August 25, 2006 was dedicated to an intercomparison of the in situ 

525 measurements on the NASA DC-8 and the British BAel46. The DC-8 made a nearly 

526 straight and level flight through a dust cloud near 2 km. Figure 18 (a) and (b) show the 

527 DC-8 altitudinal flight tracks, and the calculated S a ratios (and the coarse number 

528 concentration) for this flight, respectively. The dust layer is quite tenuous with maximum 

529 coarse number concentrations ~ 20 cm' 3 . Scattering coefficients varied from 50 to 75 

530 Mm' 1 on intercomparison legs near 19 deg N latitude. Condensation Nuclei (CN) 

531 concentrations in the dust layers were fairly low, around 300 cm' 3 , while CO mixing 

532 ratios were ~85 ppbv and RH was -50 to 60%. Both the CCN and CO concentrations 

533 infer the predominance of dust in the aerosol layer. The profiles of the 532- and 1064-nm 

534 S a ratios shown in Figure 18 for this dust layer are quite consistent with means of 38.0 ± 

535 2.5 sr, and 48.7 ± 3.2 sr, respectively. The small standard deviations in both S a ratios 

536 indicate that the layer remains very uniform with respect to this optical property. 

537 

538 Size distribution measurements were made during the NAMMA DC-8 Flight 8 on August 

539 26, 2006 of a low density dust layer between 0.6 and 1.6 km. Profiles of S a ratios 

540 calculated from these measurements are shown in Figure 19. The calculated values are 

541 39.0 ± 1.5 sr and 45.9 ± 2.2 sr at 532-, and 1064-nm, respectively. 
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542 This layer is optically thinner than the dust layer encountered on August 19, 2006. 

543 Nevertheless, the S a ratios determined by T-Matrix calculation for this layer and the direct 

544 measurement for the denser layer on August 19, 2006 (40.1 and 50.9 sr at 532 and 1064 

545 nm, respectively) are quite close. The calculated S a ratio at 532 nm is also consistent with 

546 the value (38.2 sr) retrieved from the CALIPSO measurements on the same day. 

547 

548 During the return flight on August 26, 2006, the DC-8 performed a stair-step descent 

549 flight consisting of two level sections and two descent sections. The first straight and 

550 level section (Leg 1 in Figure 20a.) was flown at an altitude of -2.25 km in a dense dust 

551 plume with coarse number concentrations near 35 particles cm" 3 . The mean S a ratios are 

552 42.4 ± 1.3 and 53.3 ± 2.0 sr at 532 nm and 1064 nm, respectively. The plume is fairly 

553 consistent as shown by the small standard deviations in the S a ratios at both wavelengths. 

554 Fig. 20(b) are profile plots of the two descent legs (leg 2 and 4 shown in plot (a) the flight 

555 path). The break in ordinate demarcates the straight and level section (leg 3) of the flight. 

556 The 532 nm and 1064 nm S a ratios for leg 2 are 46.6 ±1.3 sr, and 52.6 ±1.7 sr, 

557 respectively. In leg 4, the aircraft has begun its descent into the marine boundary layer 

558 and there is a sharp decline in the coarse number concentration. The S a ratio at 532 nm 

559 has also dropped considerably signifying a change in the aerosol composition. The 

560 calculated S a ratios for leg 4 are 32.8 ±1.5 sr, and 51.4 ± 10.8 sr, at 532 and 1064 nm, 

561 respectively. Note that for leg 4, the 532 nm S a ratios are fairly consistent while the 1064 

562 nm values are quite noisy and decrease at lower altitudes. The decreasing trend in lidar 

563 ratios at 1064 nm with altitude corresponds to a decrease in the aerosol coarse mode 

564 number concentration. 
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565 


Figure 21 shows the flight path for the in-situ measurements on August 20, 2006. The 

566 time is in seconds after midnight UTC. Plot (b) shows a profile plot of the 532 nm and 

567 1064 nm S a ratios observed during the descent phase through a dust layer extending from 

568 an altitude of 1.6 km to 4.8 km. The 532 nm and 1064 nm S a ratios for this dust layer are 

569 40.8 ± 3.2 sr, and 51.6 ± 3.8 sr, respectively. Plot (c) shows the profiles of the 532 nm 

570 and 1064 nm S a ratios observed during the ascent phase through the dust layer. The 532- 

571 nm and 1064-nm S a ratios for this dust layer are 42.8 ± 3.0 sr, and 51.8 ± 3.4 sr, 

572 respectively. The similarity of the vertical extent and the S a ratios at each wavelength 

573 suggests that the same dust layer was sampled during both the ascent and descent legs of 

574 the flight. As noted before, these layers have spatially and temporally uniform S a ratios 

575 and perhaps by inference, fairly constant compositions, and are geometrically quite 

576 stable. 

577 

578 Figure 22 is a histogram of all the 532-nm and 1064-nm S a ratios (-1100 points) 

579 determined using the size distributions measured during NAMMA for this study. There is 

580 very little overlap of the two nearly normally distributed S a ratios. The 532-and 1064-nm 

581 mean S a ratios (± one standard deviation) are 39.1 ± 3.5 sr and 50.0 ± 4.0 sr, respectively. 

582 The 532-nm values ranged from 30 to 53 sr and the 1064-nm values ranged from 32 to 66 

583 sr, in both cases within the estimated ranges of 10 - 110 sr for all aerosol types 

584 [. Ackermann , 1998; Anderson et al., 2000; Barnaba and Gobbi , 2004], Figure 23 is a plot 

585 of the frequency distribution of the ratio of S a ratios, i.e., S a (1064 nm) /S a (532 nm), a 

586 parameter used in the lidar ratio determination scheme outlined in Cattrall et al. (2005). 

587 The plot shows that for the Sahara dust sampled during NAMMA there is very little 
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588 spread in the ratio of S a values. The mean 1064 nm S a ratio is about 30% larger than the 

589 mean 532 nm value for this Saharan dust with a standard deviation of 10%, i.e., S a (1064 

590 nm) /S a (532 nm) = 1.3 ± 0.13. Since the S a ratio is an intensive property of the aerosol, 

591 its ratio is also an intensive property. The small spread in the ratios of S a denotes that the 

592 Saharan dust aerosol observed during this period is quite consistent at least in its size 

593 distributions. To explore the effects of varying composition (refractive indices) and shape 

594 (aspect ratios) on the T-Matrix calculations, both of which were not directly measured for 

595 this study, we use the uncertainty analysis described in the next section. 

596 Table 2 summarizes the S a ratios obtained for Saharan dust aerosols during NAMMA 

597 using the three methods. The 532-nm values are fairly consistent while there is a 

598 somewhat wider spread in the 1064-nm values. This is particularly interesting because 

599 these measurements were made on various days and at various locations. 

600 
601 
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Table 2. Summary of S a ratio measurements and calculations 



S a (532 nm) 

S a (1064 nm) 


sr 

sr 

Constrained method 



08/19/2006 

35.7 

NA 

Transmittance + 2-color methods 



08/19/2006 

38.9 

53.5 

08/19/2006 

41.4 

50.9 

08/19/2006 

38.1 

54.2 

08/19/2006 

41.0 

45.8 

08/19/2006 

41.0 

50.2 

08/26/2006 

38.2 

NA 

09/01/2006 

39.8 

56.0 

Mean 

39.8 ± 1.4 

51.8 ±3.6 

T-Matrix using NAMMA size 



distribution measurements 



08/19/2006 

34.3 ±2.0 

50.2 ±5.7 

08/25/2006 

38.0 ±2.5 

48.7 ±3.2 

08/26/2006 

39.0 ± 1.5 

45.9 ±2.2 

08/26/2006 

42.4 ± 1.3 

53.3 ±2.0 

08/26/2006 

46.6 ± 1.3 

52.6 ± 1.7 

08/26/2006 

32.8 ± 1.5 

51.4 ± 10.8 

08/20/2006 

40.8 ±3.2 

51.8 ±3.4 

08/20/2006 

42.8 ±3.0 

51.8 ±3.4 

Mean 

39.1 ±3.5 

50.0 ±4.0 

Noise Laden T- Matrix Simulation 



Mean 

39.4 ±5.9 

56.5 ± 16.5 


604 

605 
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606 

607 

608 

609 

610 
611 
612 

613 

614 


615 


616 

617 

618 

619 

620 
621 
622 

623 

624 


6. Uncertainty Analysis of the modeled S a ratios 

In this section we attempt to propagate the uncertainty in the input variables to the 
calculated S a ratio and determine the overall uncertainty in the calculated dust 532-nm 
and 1064-nm S a ratios using a generalized analytical uncertainty equation. The 
uncertainty equation is given by the Taylor series of the deviations (y -y°) of the output 
(y) from its nominal value (y°) and is expressed in terms of the deviations of the i inputs 
(x - x,° ) from their nominal values. As in Morgan and Henrion [1990], for the first three 
terms, the uncertainty is, 


y-y°=X ( x i - x i) Sy 


i=l 


dx- 


+ 


x° 


1 N N ^2 


i=l j=l 

N N N 


dx { dx. 


+ 


1 IN IN IN 


i=l j=l k=l 


+ ... 


X° 


I 

II (5) 

III 


The subscripts X° denote derivatives evaluated at the nominal values. Assuming that 
there are no covariances between the input variables, all the mixed derivatives in eq. (5) 
would equal zero. The only terms that would not be zero are the first term and the terms 
with i= j and i = j = k in summations II and III, i.e., second and third derivatives of the 
variables, respectively. Since the covariances are not known we cannot make the 
assumption that they are negligible. Given the uncertainties in the variables from which 
the S a ratio is calculated, the uncertainty in S a can be estimated without making any 
assumptions about covariances between inputs. To accomplish this we use Latin 
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625 


Hypercube Sampling [LHS, Iman and Conover, 1980], a statistical sampling method in 

626 which a distribution of plausible scenarios of parameter values is generated from a 

627 multidimensional distribution. Unlike classical Monte Carlo sampling methods, LHS 

628 precludes duplication by requiring that each square grid containing sample positions has 

629 only one sample in each row and each column. For this study, we generate 500 variables 

630 for each of the seven uncertain parameters used in the calculation of S a . We then 

63 1 randomly combine these variables to yield 500 instances or events. Each of these events 

632 has a very high probability of yielding a unique S a ratio. We then perform standard 

633 descriptive statistics on the S a values. 

634 

635 The mean and standard deviations of these values is an estimate of the nominal value and 

636 uncertainty of dust S a . We use the nominal (or central) values in Table 3, suggested by 

637 the studies referenced in section 2, to generate the random scenarios. The fine mode radii, 

638 coarse mode radii, and imaginary refractive indices are log normally distributed. The fine 

639 and coarse geometric standard deviations (GSD), the real refractive indices, and the 

640 aspect ratios are normally distributed. 

641 

642 The results obtained by using this method do not assume that the input variables are 

643 independent of each other or that S a ratio is linear in the individual input variables, i.e., 

644 the second- and higher-order derivatives in Equation 5 are not necessarily equal to zero. 

645 The statistics of the resulting S a ratios provide an uncertainty envelop of the S a ratio 

646 estimates based on the uncertainty of the inputs. Moreover, the results can be used to 

647 explore the sensitivity of the S a ratios at each wavelength to the various aerosol 
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648 properties. In Figure 24, the 532-nm S a ratios are well constrained with a standard 

649 deviation of 15% of the mean after perturbing nominal input values shown in Table 3 and 

650 similar to the distribution shown in Figure 22. The 1064-nm S a ratios are much more 

651 sensitive to the addition of noise as shown by the wide spread of 1064-nm S a ratios in 

652 Figure 24. A parameter that has a significant impact on the S a ratios is the complex 

653 refractive index. Figure 25 is a 2-D histogram of the 532- and 1064-nm S a ratios as 

654 functions of the real and imaginary parts of the refractive indices. The range of these 

655 values is from 0.00067 to 0.006 for the imaginary part and 1.45 to 1.55 for the real part. 

656 The 1064-nm S a ratios are sensitive to changes in the refractive index throughout the 

657 ranges of these variables. The 532-nm S a ratio is sensitive to changes in the complex 

658 refractive index at the lower ranges (Figures 25a and 25c). The highest density of 532-nm 

659 S a ratios are found in midranges of both the real and imaginary part, and these values are 

660 insensitive to changes in the complex refractive index. Note that this variation is a total 

661 derivative, i.e., all parameters are allowed to vary independently for the scenarios. The 

662 1064-nm S a ratios decrease with the complex index of refraction almost monotonically 

663 through the range of values (Figures 25b and 25d) 

664 

665 
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666 Table 3. Ranges of the variables used to generate 500 random combinations of inputs for 

667 T-Matrix calculations. The central values of the size distributions are based on the mean 

668 values of the dust layer observed on 08/19/2006 during the NAMMA campaign 

669 


Parameter 

Minimum 

Maximum 

Geometric Fine Radius 

0.0216 

0.194 

Geometric Coarse Radius 

0.188 

1.69 

Fine GSD 

1.60 

1.80 

Coarse GSD 

1.50 

1.70 

Real Refractive Index 

1.45 

1.55 

Imaginary Refractive Index 

0.00067 

0.006 

Aspect Ratio 

1.70 

2.30 


670 

671 7. Conclusion 

672 We have determined the S a ratios of Saharan dust layers using three methods: 

673 transmittance constraint technique for lofted layers, in-situ extinction profile constraint 

674 method, and T-Matrix calculations. We found quite robust S a ratios at 532 nm and a 

675 slightly wider spread at 1064 nm. The three methods yielded 532 nm and 1064 nm S a 

676 ratios that are quite close. S a ratios of 39.8 ± 1.4 sr and 51.8 ± 3.6 sr at 532 nm and 1064 

677 nm, respectively, were produced using the transmittance and 2-color methods applied to 

678 CALIPSO measurements of Saharan dust lofted layers. T-Matrix calculations applied to 

679 size distributions measured aboard the NASA DC-8 during NAMMA yielded S a ratios of 

680 39.1 ± 3.5 sr and 50.0 ± 4 sr at 532 nm and 1064 nm, respectively. The measured 

681 extinction profile obtained by the aggregate of nephelometer measurements of total 
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682 scattering coefficient and the PSAP measurements of the absorption coefficient was used 

683 to constrain the inversion of the CALIPSO measurements. This technique yielded a 532- 

684 nm S a ratio of 35.7 sr for a dust layer and 25 sr for the aerosol in the marine boundary 

685 layer. We perturbed seven microphysical and chemical properties of the dust aerosol and 

686 computed the S a ratios by randomly combining distributions of these parameters to assess 

687 the uncertainty in the S a calculation. This uncertainty simulation generated a mean (± 

688 uncertainty) S a of 39.4 (± 5.9) sr and 56.5 (± 16.5) sr at 532 nm and 1064 nm, 

689 respectively, corresponding to percent uncertainties of 15% and 29%. The ratio of the S a 

690 ratios, S a (1064 nm) /S a (532 nm) = 1.3 ± 0.13, and nearly normally distributed. The 

691 simulation revealed that S a is insensitive to the middle range of refractive indices at 532 

692 nm and nearly monotonically dependent on the refractive indices at 1064 nm. This 

693 explains the observed robust S a ratios at 532 nm and relatively wide spread of S a at 1064 

694 nm 

695 This study examined a wide range of dust loadings in various locations within -1250 km 

696 of Saharan Desert source regions on different days. The S a ratios do not change 

697 significantly with dust loading, suggesting the dust microphysical and chemical 

698 properties do not vary appreciably. The profile of S a ratios at 532 nm are very close to 40 

699 sr using different methods and on different days. There is a wider variation in the 1064 

700 nm S a ratios for all methods. It is possible that the variation of S a at 1064 nm is driven by 

701 a higher sensitivity to changes in the complex index of refraction at the 1064 nm 

702 wavelength. The ratio of S a ratios at the two wavelengths [S a (1064 nm)/S a (532 nm)] is 

703 1.3 ±0.13. 
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987 Figure 1. (a) A Two-dimensional probability plot of NAMMA in situ size distributions 

988 measured in a dust layer over a 22 minute period (-260 size distributions) during Flight 4 

989 on August 19, 2006. The mean measured size distribution is denoted by the yellow line, 

990 and (b) the best approximate bimodal log-normal particle size distribution (blue squares) 

991 with fine (coarse) mean radius of 0.0648 (0.5627) pm and a fine (coarse) geometric 

992 standard deviation of 1.696 (1.572) fitted to the average measured size distribution (red 

993 squares). 
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995 



996 Geometric Standard Deviation Geometric Standard Deviation 

997 Figure 2. The distributions of the fine (a) and coarse (b) geometric mean radii and the 

998 fine (c) and coarse (d) geometric standard deviations for a dense elevated dust layer 

999 measured during a flight on August 25, 2006. 
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Figure 3. The (a) 532-nm and (b) 1064 nm Saharan dust phase functions calculated from 
the size distributions measured aboard NASA DC-8 during Flight 4 on August 19, 2006. 
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1006 



1007 Seconds after 00:00 UT 

1008 Figure 4. The DC-8 altitudinal flight path showing the segment during which the in situ 

1009 measurements were made. The direct CALIPSO underflights are the five segments (a, b, 

1010 c, d, and e) at the highest altitudes. The time is in seconds after midnight UTC. 
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Figure 5. CALIPSO flight tracks superimposed on (a) the MODIS optical depth and (b) 
the MOPPITT CO concentrations. The DC-8 flight track for August 19, 2006 is shown 
by the black lines in the middle of both figures. The blue dots in (a) denote the portion of 
the DC-8 flight segment during which the in-situ measurements were taken. Daytime 
CALIPSO tracks run from the southeast to the northwest while the nighttime tracks run 
from northeast to the southwest. 
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1021 



1023 Figure 6. Browse image of the CALIPSO flight track on August 19, 2006 shown in blue 

1024 in the inset map and the location of the near coincidence of the CALIPSO flight track and 

1025 the DC-8 flight 4 during NAMMA is denoted by the two white lines near 15N, 19W. In 

1026 the image, colors above the surface are coded as follows: blue indicates clear sky, white 

1027 indicates clouds, and warm colors (green, yellow, orange, and red) indicate aerosols. The 

1028 white parallel lines bound the CALIPSO profiles closest to the in-situ measurements. 
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1032 Figure 7. A comparison of the profiles of the in situ extinction coefficient derived during 

1033 the DC-8 ascending segment as shown in Figure 4 and the extinction coefficient retrieved 

1034 from the CALIPSO measurement averaged over the region indicated by the two white 

1035 lines in Figure 6. The S a ratios were determined using the in situ extinction as an 

1036 additional constraint. 
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1039 Figure 8. Aerosol scattering ratio (p a / p m at 815 nm) image measured by the Lidar 

1040 Atmospheric Sensing Experiment (LASE) aboard DC-8 during the ascending segment 

1041 passing through a Sahara Air Layer (SAL) and the underflight of CALIPSO as shown in 

1042 Figure 4. The DC-8 point of nearest approach to the CALIPSO flight track is near 15°N- 

1043 19°W. 

1044 


49 


1045 

1046 

1047 

1048 

1049 

1050 

1051 


10.0- 
8 . 0 - 
'E 6 . 0 - 

Q) 

-d 4.0- 

ZJ 

4 ^ 

< 2 . 0 - 

0 . 0 - 

lO.Ch- 

8 . 0 - 

^ 6.0- 
& 

"B 4.0- 

+3 

-f- 

^ 2 . 0 - 

O.CLz 

Lat 

Lon 




-8.8 -6.2 -3.7 -1.1 1.5 4.0 6.6 9.2 11.7 14.3 16.8 

-13.7 -14.3 -14.8 -15.4 -15.9 -16.5 -17.0 -17.6 -18.1 -18.7 -19.3 


Marine Dust 


Polluted 
Continental 


„. Polluted 

Clean hll<+ Smoke 

Continental 


Figure 9. Products of the CALIPSO (a) cloud-aerosol discrimination and (b) aerosol 
classification schemes for the data in the browse image shown in Figure 6. The cloud 
aerosol discrimination shows the extensive aerosol layer at 3 km and the aerosol 
classification image shows that this aerosol is a Saharan dust layer to the north and a 
polluted dust and smoke layer to the south. 
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Figure 10. (a) Browse image of the orbit segment denoted by the blue track in the inset 
map and (b) an expansion of the dust layer region of the plot marked by yellow dots in 
(a). In both cases the color code is as in Figure 6. The red dashed oval shows the elevated 
(layer top ~ 5km) dust layer. Aerosol optical depths at 532 nm for all five segments are 
shown on the figure and are near 0.3. 
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Figure 11. CALIPSO flight tracks superimposed on the (a) MODIS optical depth and the 
(b) MOPPITT CO concentrations. The CALIPSO flight tracks are colored cyan and grey 
in the MODIS and MOPPITT images, respectively. The DC-8 flight track for August 26, 
2006 is shown by the heart shaped drawing in the middle of both figures. The CALIPSO 
underflight segment is the portion of the flight track in alignment with the CALIPSO 
flight track in cyan or grey in (a) or (b), respectively. The DC-8 in situ data for this 
obtained during the descent flight segment shown by a red dashed tilted oval in (a). 
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Figure 12. Browse image of the CALIPSO flight track shown in green in the inset map 
and the location of the near coincidence (in space) of the CALIPSO flight track and the 
DC-8 flight 8 on August 26, 2006. The color code is as in Fig. 6. The white lines denote 
the CALIPSO data used for comparison with the DC-8 measurements. 
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1078 Figure 13. A comparison of the NAMMA in situ extinction profile and the extinction 

1079 profile calculated from CALIPSO’s attenuated backscatter measurement using 532 nm S a 

1080 ratio of 38.2 sr obtained using the transmittance method. The two measurements both 

1081 taken on August 26, 2006 are offset by 2 hours and ~1250 km. 
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Figure 14. CALIPSO flight tracks superimposed on the (a) MODIS optical depth and the 
(b) MOPPITT CO concentrations. The CALIPSO flight tracks are colored cyan and grey 
in the MODIS and MOPPITT images. The DC-8 flight track for September 1, 2006 is 
shown by the black lines in the middle of both figures. 
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Figure 15. Browse image of the CALIPSO flight track shown in cyan in the inset map 
and the location of the near coincidence (in space) of the CALIPSO and the DC-8 flight 
tracks on September 1 during NAMMA. The color code is as in Fig. 6. The white lines 
subtend the region averaged for the extinction comparison. 
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1099 Figure 16. The NAMMA/CALIPSO extinction comparison for Flight 10 on September 1, 

1 100 2006. Note that the centroid of the layer has shifted upwards by about 1.5 km in between 

1 101 the CALIPSO measurement and NAMMA sampling. The CALIPSO measurement is ~15 

1 102 hours earlier and nearly 750 km away. 
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Figure 17. A profile of the 1064 and 532-nm S a ratios and integrated number distribution 
of coarse particles (0.583 um - 6.264 um) taken from a 20-minute segment of the APS 
measurements during Flight 4 on August 19, 2006. The dust base and top are shown by 
the lines at altitudes of 2.5 km and 4.6 km, respectively. The mean 532-nm and 1064-nm 
S a ratios for the layer is 34.3 ± 2.0 sr and 50.2 ± 5.7 sr, respectively. 
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1112 Figure 18. (a) The flight track of the DC-8 on August 25, 2006. The straight and nearly 

1113 level portion is through a dust layer at ~2 km. (b) The calculated 532- and 1064-nm S a 

1114 ratios and coarse number particle concentration for the straight and level portion of this 

1115 flight. 
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Figure 19. A profile of the 1064- and 532 -nm S a ratios, and integrated number 
distribution of coarse particles during Flight 8 on Aug. 26, 2006. The mean 532-nm and 
1064-nm S a ratios for the layer is 39.0 ± 1.5 sr and 45.9 ± 2.2 sr, respectively. 
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1129 Figure 20 . (a) Flight path and S a ratios and the particle number concentration for the 

1130 return flight on August 26, 2006. The figures are plotted by temporal profiles for the 

1131 straight and level legs, and vertical profiles for descent legs, (b) is temporal profile of 

1132 straight and level leg 1, (c) vertical profiles of descent legs 2 and 4, and (d) temporal 

1 133 profile of straight and level leg 3. 
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1136 Figure 21. The (a) flight path and S a ratios and the particle number concentration for the 

1 137 flight on August 20, 2006 plotted vertical profiles of (b) the descent leg and (c) the ascent 

1138 leg. 
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Figure 22. An ensemble of all the S a ratio estimates (-1100 records) based on the 
NAMMA in-situ size distributions of dust layers and T-Matrix calculations. The average 
S a ratios are 39.1 ± 3.5 sr and 50.0 ± 4 sr at 532 nm and 1064 nm, respectively 
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Figure 23. The ratio of the S a ratios [S a (1064 nm)/S a (532 nm)] derived from the size 
distribution data used for this study. 
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Figure 24. The distribution of the S a ratio calculated from a Latin Hypercube Sample 
(LHS) of 500 events bounded by the ranges shown in Table 3. The normal means (± 1 
standard deviation) of the 532- and 1064-nm S a ratios are 39.4 ± 5.9 sr and 56.5 ± 16.5 sr, 
respectively. 


67 




1160 

1161 

1162 

1163 

1164 

1165 



Imaginary Part of the Refractive Index 


90 


70 

60 

50 

40 

30 











(b) 

■l 

1 




















V 















I 





1 










% 

- 











V 

1 L 
















0.001 0.002 0.003 0.004 0.005 

Imaginary Part of the Refractive Index 



7 

6 

5 

4 

3 

2 




1.46 1.47 1.48 1.49 1.5 1.51 1.52 1.53 1.54 

Real Part of the Refractive Index 



1.46 1.47 1.48 1.49 1.5 1.51 1.52 1.53 1.54 

Real Part of the Refractive Index 


Figure 25. Two-dimensional frequency plots showing the variation of the lidar ratios with 
real and imaginary parts of the complex refractive indices for 500 randomly generated 
events discussed in section 6. 
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